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The size of the vortex core in a clean superconductor is strongly temperature dependent and 
shrinks with decreasing temperature, decreasing to zero for T ^ 0. We study this so-called Kramer- 
Pesch effect both for a single gap superconductor and for the case of a two gap superconductor 
using parameters appropriate for Magnesium Diboride. Usually, the Kramer-Pesch effect is absent 
in the dirty limit. Here, we show that the Kramer-Pesch effect exists in both bands of a two gap 
superconductor even if only one of the two bands is in the clean limit and the other band in the 
dirty limit, a case appropriate for MgB2. In this case an induced Kramer-Pesch effect appears in the 
dirty band. Besides numerical results we also present an analytical model for the spatial variation 
of the pairing potential in the vicinity of the vortex center that allows a simple calculation of the 
vortex core radius even in the limit T ^ 0. 

PACS numbers: 74.20.-z, 74.25. Bt, 74.70.Ad 



I. INTRODUCTION 

In 1974 Kramer and Pesch theoretically discovered 
an unusual core shrinkage of an isolated vortex for de- 
creasing temperature which they did not only observe in 
numerical solutions of Eilenberger's equations but also 
proved analytically [j. While all other lengthscales de- 
scribing the superconducting state, especially the Lon- 
don penetration depth and the coherence length, reach a 
saturation value with decreasing temperature, this inves- 
tigation introduced a new lengthscale in the discussion of 
the vortex state of clean superconductors. It can be de- 
fined as the inverse of the slope of the pairing potential at 
the vortex center and describes not only the spatial vari- 
ation of the gap function but also the maximum height 
of the supercurrent density around the vortex center and 
thus measures the size of the vortex core. In a clean su- 
perconductor without impurity scattering this length ^„ 
decreases linearly with temperature tending to zero for 
T while in superconductors with larger scattering 
rates it reaches a saturation value. A detailed study of 
the impurity effect on the vortex core shrinkage can be 
found in Ref. 0. 

It appears well established both experimentally and 
theoretically that the new superconductor MgB2 is a 
two gap superconductor, possessing two superconduct- 
ing gaps of significantly different size on disconnected 
parts of the Fermi surface There exist other 

I compounds, in which two gap superconductivity is be- 
lieved to be reaHzed, however, at present MgB2 is the 
clearest example. In the case of such a multi band su- 
perconductor with several distinct gap values there exist 
different lengthscales that describe the spatial variation 
of the quasiparticle excitations corresponding to the dif- 
ferent gap values. Therefore, we also expect different 
lengthscales for the description of the increase of the gap 
functions near the center of an isolated vortex in each 
band. It is obvious that these lengthscales are not inde- 



pendent of each other if there is a coupling of the different 
bands via the pairing interaction. In this work we want 
to examine numerically the vortex core shrinkage for a 
two gap superconductor at the example of Magnesium 
Diboride. We also present an analytical gap model that 
well describes the low temperature behaviour of the core 
size ^t,. In the last section we want to discuss a model 
of a two gap superconductor assuming large and small 
scattering rates in the different bands and its infiuence 
on the size of the vortex core. 



II. THE SINGLE BAND CASE 

In this section, we study the vortex core structure of a 
standard (single band) superconductor. In the first sub- 
section IlII All , we describe how the Riccati parametriza- 
tion formahsm of the Eilenberger theory can be used 
to numerically calculate the structure of an isolated 
Abrikosov vortex in the clean limit. In the second subsec- 
tion dllBp , we present an analytical model for the pairing 
potential in the vicinity of an isolated Abrikosov vortex 
and explain how it can be used to calculate the tem- 
perature dependence of the vortex core size. The third 
subsection ijllCII appHes to the dirty limit. We point out 
how the theoretical description simplifies in this special 
case and discuss the results for the vortex core size. 



A. Numerical Calculations 

In order to calculate the pairing potential A(r, T) for 
a certain region in real space, one has to find a solution 
to the gap equation: 

A(f,T) = VN{0)27rT ^ (/(r, fcf, i£n))f s 

0<£„<Wc 

(1) 
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According to Refs. 0, the quasiclassical propagator 
f{r,kp,ien) can be expressed in terms of two complex 
quantities a{x) and h(x), called the Riccati amplitudes: 



f{f{x), kp^iSn) 



2a{x) 



1 + a{x) b{x) 



(2) 



The Riccati amplitudes a{x) and b{x) are in turn solu- 
tions to the Riccati differential equations 

hvpd^ a{x) + [2e„ + A''{f{x)) a{x)] a{x) - A{r{x)) = 

(3a) 

hvpd^ b{x) - [2e„ + A(f(a;)) b{x)] b{x) + A^(f(a;)) = 

(3b) 

The Riccati differential equations have to be solved along 
real space trajectories r{x) pointing in the direction of the 
Fermi wave vector vpikp) using the modified Matsubara 
frequencies ?£„ — iSn + (e/c) vp ■ A{r{x)). The quasiclas- 
sical Green's function f{f,kF,i£n) calculated with the 
Riccati amplitudes according to equation IpJ then has to 
be averaged over the Fermi surface of the superconduct- 
ing material (denoted with (• • • )fs)- 

As initial values for the integration of the Riccati differ- 
ential equations the bulk values of the Riccati amplitudes 
have to be used (e„ > 0): 



a(— oo) 
6(-l-oo) 



A(-cx)) 



e;i + |A(~oo)r 
At(+oo) 



(4a) 
(4b) 



The quasiclassical Green's function /(r, kF,iSn) is, via 
the Riccati amplitudes, a function of the pairing poten- 
tial. A straightforward way to find a solution to the 
self-consistency problem posed by the gap equation is 
to choose a certain initial configuration for the pairing 
potential A(r, T) and then improve this configuration it- 
eratively. 

Since we will apply the methods described in this chap- 
ter to study the temperature dependence of the vortex 
core size of an isolated Abrikosov vortex, it is convenient 
to define, at this point, a characteristic lengthscale for 
the size of the vortex core. We take up the proposal for 
a definition by Hayashi et al. 0: 



5A(r) 



dr 



A{r = oo, T) 



(5) 



This characteristic lengthscale ^i, corresponds to the in- 
verse slope of the pairing potential A(f, T) at the vortex 
center r = 0, normalized to the bulk value of the pairing 
potential A{r — oo, T) at the corresponding temperature 
T and has to be distinguished from the bulk coherence 
length ^oo = hvp/Air = oo, T). 

Once a self-consistent solution for the pairing potential 
is found, the local quasiparticle spectrum N{r, E) (local 




Figure 1: The quasiparticle spectrum in the single band case, 
plotted for several distances from the vortex core at a temper- 
ature of T = 0.5 Tc. The distances span from r = (rearmost 
curve) to r = 5 ^oo (foremost curve) in steps of 0.2 ^oo. For 
these plots, an imaginary part (broadening) of 5 = 0.1 was 
used. 



density of states) can be calculated in the following man- 
ner: 



N{r, E) 





'1 


— ab 


{Re 




1 


+ ab 



(6) 



FS 



In Fig. m we show numerical results for the quasipar- 
ticle spectrum in the single band case. The curves were 
obtained evaluating Eq. lO based on self-consistently cal- 
culated data for the pairing potential A(r, T) and a cylin- 
drical Fermi surface has been assumed. A distinct zero 
energy peak represents Andreev bound states in the vor- 
tex core. For increasing distance from the vortex core, 
this peak splits and finally merges into the gap edges. 
Similar results have been obtained before. Gygi and 
Schliiter solved the Bogoliubov-de Gennes equations as 
an eigenvalue problem to study the electronic structure 
of a vortex line, see Ref. @. Ichioka et al. used the qua- 
siclassical Eilenberger theory to examine the same prob- 
lem, see Ref. 

Although the Riccati parametrization of the Eilen- 
berger propagator offers a convenient and stable method 
to solve the Eilenberger equations and to calculate the lo- 
cal density of states, the effort of numerical calculations 
increases under certain circumstances. Firstly, as the 
lengthscale on which the pairing potential A(f, T) varies 
is getting shorter, the number of grid points on which 
the gap equation has to be solved in real space grows. 
Secondly, the number of terms in the sum over Matsub- 
ara frequencies J2o<6 <uj grows linearly with decreasing 
temperature. Thirdly, the averaging over all Fermi wave 
vectors vp that appear on the Fermi surface becomes la- 
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borious for complicated Fermi surface topologies. 

Taking all these points into account, it is obvious that 
an analytical approach should be considered whenever 
possible, in particular if one wants to study the T ^ 
limit. In the next chapter, this will be done for an isolated 
Abrikosov vortex. 



B. Exact Solution for an Analytical Gap Model 



In the innermost part of an Abrikosov vortex, the mod- 
ulus of the pairing potential grows linearly with increas- 
ing distance from the vortex center. As it approaches 
the bulk value, it saturates with a certain profile. Su- 
perimposed to the behaviour of the modulus of the pair- 
ing potential, the phase varies by 27r around the vortex. 
This leads to the general form A(r) ~ re"^ for the in- 
nermost part of the vortex, while in the outermost part, 
the behaviour A(r) ~ e*"^ holds (r and 4> are polar coor- 
dinates in real space). Thus, the following simple analyt- 
ical model for the pairing potential in the vicinity of an 
isolated Abrikosov vortex is useful: 



A(f,T) 



Aoo(T) 

Aj(T)e^"* 



,r <£.v 



(7) 



Here, Aoo(T) denotes the value of the pairing potential 
in the bulk at a given temperature T; is the length- 
scale on which the pairing potential rises from zero in 
the center of the vortex to the bulk value Aoo{T) and 
coincides with the previous definition ©. According to 
Ref. 0, every real space trajectory f{x) is parametrized 
by an impact parameter y and the position along the tra- 
jectory X. In our special case, the impact parameter y is 
given by the distance from the trajectory to the vortex 
core. The position along the trajectory x is measured 
relative to the point closest to the vortex core. It follows 
from what has been said that 



pi4> 



{x + iy) e 



(8) 



along a chosen trajectory specified by 6*, the angle cor- 
responding to the direction of the Fermi velocity in the 
a6-plane. Within the scope of this model, the pairing po- 
tential consists of a vortex core with a linear profile and 
a phase vortex in the circumference. 
In the next two subsections, we will show how an analyti- 
cal solution can be obtained for the quasiclassical Green's 
function in the special case of this particular pairing po- 
tential profile. We will then use this analytical solution 
in order to determine the lengthscale S,^ self-consistently 
with the gap equation. 



1. Analytical Solution for r < 

In general, the linearized Bogoliubov-de Gennes equa- 
tions, often referred to as Andreev equations, read 



ax 



u(x) 




v{x) 





i£n{x) -A(a;) 
A^(a;) -Unix) 



u{x) 
v{x) 



(9) 



According to Ref. the solutions of Q are connected 
to the Riccati amplitude a (a;) via 



a{x) 



u{x) 
v{x) 



(10) 



whereas the Riccati amplitude h{x) can be obtained ap- 
plying the symmetry relation 



b{x) 



-a{—x) e 



-2ie 



(11) 



which is valid if only a single vortex line is present. 
It should be noted that we do not solve the linearized 
Bogoliubov-de Gennes equations ijHl as an eigenvalue 
problem, but as an inital value problem, see Ref. 
For the linear profile, Eq. Q reads 



e'2 ' 



i£nh - {xT2 + yfi) e *2^3 1 



where n are the Pauli matrices. Applying the gauge 
transformation 



*i(x) 
*2(x) 



u{x) 
v{x) 



(13) 



and introducing the spinor $(x) via 

A^ 



hvpdx - En fa 



{xf2 + yfi) 
oe^f^i o ^x) (14) 



one obtains the following two decoupled differential equa- 
tions for the components $(a;) = ($_(a;), $+(a:))"^ 



ScaHng the variables x 

acteristic lengthscale ^ 
introducing 



o^^{x) = Q (15) 

I to the char- 
- Ji^ ^ith Coo = and 



and y 



$±(a;) = 



hvp 

4 U Ao< 



2 



(16) 



(17) 
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Eqs. 1115(1 assume the following compact and dimension- 
less form: 



dl 



1 -2 



$±(x) 







(18) 



The general solution of this second order differential 
equation is given by a linear combination of the parabolic 
cylinder functions U{c±,x) and V[c±,x) fTll|. We choose 
the following ansatz 

$_(s) = A_y(c_,-x)e(-5) + B_y(c_,x)e(x) 

l'+(x) = A+C/(c+,-i)e(-5:) +B+f7(c+,x)e(i) 

(19) 

1 for a; > 
for S < 



with 6(x) 



in order to construct the Riccati amplitudes a (a;) and 
h{x). The connection between the Riccati amplitude a(x) 
and the solutions of the Andreev equation II12II , given by 
ltTn)l . now reads 



a{x) 



*2(a;) 



(20) 



Let us briefly comment on the construction of the Ric- 
cati amplitudes a{x) and h(x), respectively. A correct 
solution firstly has to solve the Riccati equations (El and 
secondly has to possess the correct bulk asymptotics I^J- 
The former is provided by construction |g| and we will 
see later on that the asymptotic behaviour obtained with 
this ansatz is indeed correct. 

Additionally, it should be noted that the ansatz l|19|l is in 
fact unique. An expansion of the solution derived from 
the most general ansatz (containing both parabolic cylin- 
der functions U(c±,x) and V{c±,x) for both functions 
^± and all x) exhibits the correct asymptotic behaviour 
if and only if the terms not contained in l(19|l vanish. 
Thus, the above form is being reproduced. 



Using the ansatz l|19|l . undoing the gauge transforma- 
tion Ill8|l and the substitutions l(14|l . Ijlfill . then applying 
some basic properties of the parabolic cylinder functions 
[m and making use of l|2()ll . one obtains the following 
expression for the Riccati amplitude a(x) for x < 0: 



a(x) 



i0 


[ (7(c_,-2) 2 \ i '■y) 


V(c_,-S)l 








U(c+,-x) 
U(c^-x) 






y(c_ -2)1 

(7(c_,-£) 








U{c+,~x) 
U(c^-x) 



(21) 



A similar expression can be obtained for a; > 0, but 
this will not be necessary for our further calculations. 
The Riccati amplitude b{x) for a; > can be constructed 
from equation l|21ll using the symmetry relation llllll . We 
now have to determine the coeflicient {^-) in order to 
complete the solution. This will be achieved by matching 
Eq. ll2T|l to the solution for r >£,v 



2. Analytical Solution for r > 

The outer region of our model for the pairing potential 
is nothing but a pure phase vortex. Unfortunately, 
there is no exact solution to the Riccati equations (pj 
for such a phase vortex for all energies and all impact 
parameters. On the one hand, an exact solution for all 
values of the impact parameter can be obtained for the 
special case iSn ^ E = |Aoo| 3. On the other hand, an 
asymptotic solution for all energies is obvious for vanish- 
ing impact parameter (In this special case, the pairing 
potential is constant along any single trajectory except 
for the phase step in the vortex core) . Since we will solve 
the gap equation for a real space position Xq that is very 
close to the vortex center (in fact we will perform the 



limit xq — > 0), the impact parameter will be vanishingly 
small in our calculations. Thus, we choose the asymp- 
totic solution that is exact for all energies and vanishing 
impact parameter. It is of the following form: 



a{x) 



£n + ^Jel + I Ac 



(22) 



Now the solution for the outer part of the vortex being 
constructed, we can demand that the solution be contin- 
uous at ^t,. This means that Eq. lf2T|l evaluated at a: = ^„ 
has to be equal to Eq. lf22|l which determines the coeffi- 
cient i,^-)- Hence, we found a solution for the Riccati 
amplitudes for all positions x and all energies which is 
exact for vanishing impact parameter y. 

We still have to check if the asymptotic behaviour of 
our solution is consistent with the bulk asymptotics of 
the Riccati amplitudes. These have been discussed in sec- 
tion Ijll A|l and are given in Eqs. Q. On the one hand, we 
insert the linear profile of the pairing potential into Eqs. 
Q and thus obtain the correct asymptotic behaviour. 
On the other hand, we move the matching point of our 
analytic solution to infinity and then expand Eq. lf2T|l in 
terms of a; ^ ±oo. Indeed, this reproduces the correct 
bulk asymptotics. With a finite matching point, par- 
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ticularly the matching point ^v, the correct asymptotic 
behaviour is provided by construction of the solution. 

Altogether, we can therefore state that the solution 
constructed with ansatz indeed exhibits the correct 
asymptotic behaviour in the bulk. 



3. Solving the Gap Equation 

In order to determine the lengthscale £,v on which the 
pairing potential rises from zero to the bulk value, we 
have to solve the gap equation 1^ self-consistently for a 
real space position xq which is close to the vortex center. 
Therefore, we have to calculate the quasiclassical Green's 
function /(r, fc^, according to Eq. |(5J| using the so- 
lutions for the Riccati amplitudes a{x) and b{x) we found 
in subsections ijllB 111 and ijllB 2|l . Since we are only in- 
terested in the lengthscale which corresponds to the 
inverse slope of the pairing potential in the vortex core 
(see Eq. @), we differentiate the gap equation (QJ at the 
vortex center xq = with respect to the position xq: 



Aoo(T) 



VN{0)27rT J2 {d.ofixo^O))FS 



(23) 



We now see that we only need to calculate the derivative 
of the Green's function /(r, kp, i£n) at the vortex center: 



dxofi^o = 0) 



d 2a 



dxa 1 + ab 



(24) 



xo=0 



We therefore only need to know the quantities a{xo — 0) 
and dxoa{xo = 0) which correspond to a vanishing impact 
parameter y. The Riccati amplitude b{x) can then be cal- 
culated using the symmetry relation . The expression 
obtained this way then has to be averaged over the Fermi 
surface (denoted by (• ■ • )fs) and inserted into the gap 
equation. The lengthscale £,viT) is, for the full temper- 
ature range, determined by this equation. Note that S,v 
still has to be determined self-consistently, because the 
Riccati amplitude a{x) depends on via Eq. Ij21|l . 

Fig. [2 shows the results from these calculations com- 
pared to those from a fully numerical self-consistent so- 
lution of the gap equation described in section IlII A|l . In 
both cylindrical Fermi surface has been used. 

The model for the pairing potential JJjl yields a qualita- 
tively correct behaviour of (T) over the full tempera- 
ture range. The lengthscale ^„ diverges for T ^ Tc and 
vanishes for T — > 0. In Fig.Elwe show the numerically de- 
termined gap profile together with the approximate one 
using the model given in Eq. Q. At low temperatures 
the slope at the vortex center is reproduced very well by 
the approximate model, while at higher temperatures de- 
viations occur. This is refiected in the behaviour of the 
lengthscale ^„ in Fig. |2l 

Based on the analytical expressions obtained in this 
chapter, it is possible to calculate in the limit T ^ 




Figure 2: The lengthscale (,v{T) obtained from a linear profile 
of the pairing potential (solid line) and from a numerical self- 
consistent solution of the gap equation (diamonds with dashed 
guideline). Additionally, the plot shows the inclination of 
in the limit T — > (dotted line). The inset is a blow-up of 
the low temperature region. 




Figure 3: The modulus of the pairing potential as a func- 
tion of the distance from the vortex center for three different 
temperatures. For each of the three different temperatures, 
we show the numerically determined self-consistent solution 
(solid line) and the linear profile with a slope calculated eval- 
uating the gap equation at a position xo ^ (dashed line). 
The matching of the slope at r = becomes very good at low 
temperatures. 



as long as the cut-off frequency ujc is finite. Therefore, 
one has to insert the analytical result for the Riccati am- 
plitude a{x) given in Eq. Ij21|l into Eq. Ij24ll making use 
of the symmetry relation ; this quantity can then be 
inserted into the differentiated gap equation ll23)l and the 
limit T — > can be performed. The result is the follow- 
ing: 



4^o 



7rFiV(0) Ac 



T 
TV 



for T ^ 



(25) 



This leads to £,v oc £,ooT/Tc and we can thus, within the 
Riccati parametrization formalism of the Eilenberger the- 
ory, reproduce the result found earlier by Kramer and 
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Pesch Pi as a limiting case of our calculations. 

It should be noted that the linear behaviour for T ^ 
is valid only for a very narrow temperature range (see the 
inset in Fig. EJ. The curve ^v{T) bends towards = 
not until very low temperatures (~ 0.05 Tc). Thus, linear 
extrapolations from higher temperatures towards T = 
should be handled with care. 



C. Comparison with Dirty Limit Calculations 

In the case of a dirty superconductor with high scatter- 
ing rates the Eilenberger equations reduce to a simple dif- 
fusionhke equation for the Usadel propagator represented 
by the Green's functions G'(r, e„) and F(r, £„) that are 
independent of the Fermi wave vector kp as was shown 
by Usadel 1121 . Following the notation of Koshelev and 
Golubov in this equation can be written as 



V 

SnF 

2 



G 



$0 



= AG- 



together with the self-consistency condition 



A(f) = VN{0) 2ttT 



0<e„<iJ, 



(26) 



where we have assumed that the Fermi surface is rota- 
tionally symmetric with respect to the applied magnetic 
field and we have reduced the equations to the single 
band case. The characteristic lengthscale is determined 
by the diffusion constant V and can be written as 



V 
2^ 



(27) 



Employing the normalization condition 

[G{r,en)f + F\r,en)F{r,en) = 1 

we can introduce a parametrization for the momentum 
averaged Usadel Green's functions G(r, £„) and F(r, e„) 
similar to the Riccati parametrization for the clean limit 
Eilenberger equations j3 ^| ■ 



G(r,e„) 



1 



F{r,e„ 



2a 



1 



(28) 



It should be noted that here we restricted to positive 
Matsubara frequencies £„ > 0. 

In the vicinity of an Abrikosov vortex we can simplify the 
equations by switching to cylindrical coordinates. If we 
assume a large Ginzburg-Landau parameter k ^ 1 and 
confine to an isolated vortex we can neglect the vector 
potential and write down the differential equation for a(r) 
as 



2£„a — V 



-drirdro) 
r 



2a{draf 1 a{l-o?) 



1 

,2 



Here r denotes the distance from the vortex center and 
the phase of the order parameter tp has been assumed to 
coincide with the polar angle (t> as A(r) = A(r)e"^. The 
parameter a{r) is related to the parameter of by 
a(r) — tan|. Inserting this relation into Eq. ll2iJ|l and 
normalizing all energies to ttTc and all lengths to £,oa we 
reproduce the differential equation for Q{r) from Ref. [T^ : 



1 



1 



d^0{r,en) + -dr9{r, e„) ^ cos 6l(r, e„) sin6'(r, e„) 

+A(r) cos 6'(r,e„) - £„ sin6l(r, e„) = (30) 
together with the gap equation 

A(r) = 1/A^(0) 27rT ^ sin6i(r,e„) (31) 

The boundary conditions for the differential equation 
H30II follow from general considerations: In the vortex 
center the anomalous Green's function F(r, e„) vanishes 
together with the modulus of the pairing potential lead- 
ing to 



(r = 0,£„) = 



(32) 



and far from the vortex center we expect the Green's 
functions to assume their bulk values and we obtain 



oo, £„) = arctan 



Ao 



(33) 



To obtain a self-consistent profile for the pairing potential 
A(r) the Usadel equations have to be solved iteratively 
together with the gap equation as described in section 
Ijll A|l . We used a relaxation method for a fast and sta- 
ble solution of the boundary value problem defined by 
Eqs. (jSni El and chose a cut-off for the Matsub- 
ara frequency summation in Eq. ll3T|l of uoc — 20nTc. 
Extracting the characteristic length ^„ from the inverse 
slope of the pairing potential in the vortex center we find 
- as expected from several other calculations, e.g. j3| - 
a saturation value of ^t, for low temperatures, see Fig. ^ 
i.e. in the dirty limit the Kramer-Pesch effect disappears. 

In Fig. 13 we show the quasiparticle spectrum calcu- 
lated as the real part of the analytical continuation of 
the Usadel Green's function G(r, —iE + 6) according to 
Eq. II3()|I . Due to strong impurity scattering, the zero en- 
ergy bound state in the vortex core is replaced by a fiat 
normal state spectrum. 



III. THE TWO GAP CASE 

In this section, we study the vortex core structure of a 
two gap superconductor. In the first subsection dlll A\i . 
we explain how the numerical procedure for the single 
band case from section ijll A|l can be generahzed in order 
to describe a clean two band system. In the second sub- 
we do the same for the analytical model 
TTbIi. In this sub- 



(l-a^)A(r) (29) 



section ijlll B 

for the pairing potential from section 
section, we also present results for a two gap supercon 
ductor in the dirty limit. 
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Figure 4: The lengthscale (.v{T) for the dirty limit approach. 
A saturation occurs at low temperatures. 




Figure 5: The quasiparticle spectrum in the dirty limit, plot- 
ted for several distances from the vortex core at a temperature 
of T = 0.5 Tc. The distances span from r — (rearmost curve) 
to r = lO^oo (foremost curve) in steps of 0.5^ca- 

A. Numerical Calculations 

In order to include two band behaviour into our con- 
siderations, the gap equation (QJ has to be replaced by a 
two band gap equation: 

A(")(f,r) = 27rT^A„„, (Mfs,^ (34) 

a' 0<en<ujc 

with a, a G {cr, tt} 

This two band gap equation can be derived starting from 
the fully momentum dependent multiband formulation of 
the quasiclassical (Eilenberger) theory and has been 
used successfully to describe the upper critical field in 
MgB2 We introduced two band indices (cr and tt) 

with respect to the usual nomenclature in MgB2 as well 



as a 2 X 2 coupling matrix Xaa' ■ 

The application of the two band gap equation re- 
quires the calculation of the quasiclassical Green's func- 
tion for the two bands separately, (r, A:f, for the 
a band and f-n{r,kF,i£n) for the tt band. Eqs. 
still apply, but the cr component of the pairing po- 
tential /\^"\f,T) has to be used for the calculation of 
fa{r,kF,i£n) and the tt component A('^)(r, T) for the 
calculation of f-n{r, kp,ien), respectively. 

Additionally, the Fermi surface topology of MgB2 has 
to be taken into account in the Fermi surface averages 
(• • • ) FS,<7 and (• • • ) FS,TT, respectively. Band structure cal- 
culations have been carried out to determine the Fermi 
surface structure as well as microscopic calcula- 

tions that revealed the distribution of the superconduct- 
ing gap on the Fermi surface Here, we use a simple 
model for the Fermi surface that has been suggested in 
Ref. ^ half-torus for the tt band and a distorted 
cylinder for the cr band, respectively. This model cor- 
rectly reproduces the Fermi surface topology and has 
proved to describe the upper critical field anisotropy in 
MgB2. This leads to the following integral parametriza- 
tions IHlii: 

(• ■ ■ )FS,a = J dkc J d(j) ■■■ (35) 

Here, (p G [0i27r] is the azimuthal angle within the ab 
plane, 9 & [f , ^] the polar angle of the torus, kc the c- 
axis component of the momentum, c the lattice constant 
in the c-axis direction, and k — 0.25 the ratio of the two 
radii of the torus [T^ . 

In the two band gap equation Ij34|l . the coupling be- 
tween the two bands occurs via the coupling matrix 
Xaa' which, in the weak coupling limit, can be calcu- 
lated from the ratio of the gap amplitudes for T ^ 
(A^VaS^(T ^ 0)) and for T ^ (A^^^/A^^\t ^ 
Tc)), from the ratio of the densities of states in the nor- 
mal (non-superconducting) state {N{OY'^^ /N{OY^^) and 

from the cut-off frequency ujc H^l- The ratio A^-'/A^' 
can be extracted from calculations of the two supercon- 
ducting energy gaps based on the Eliashberg formalism 
as well as the ratio 7V(0)('"V^(0)^''^ The values 

A^Va^^(T ^ 0) = 3, A^VA^^(r ^ Tc) = 4 and 
7V(0)('^V^(0)^''^ = 0.734 have been used in our calcula- 
tions. 

The general form of the two band gap equation Ij34|l leads 
to the complication that the two components of the pair- 
ing potential, A('^)(r;T) and A(^)(f;T), have to be calcu- 
lated simultaneously. 

In the case of a two band superconductor, the spectral 
distribution of quasiparticles N{f, E) is the weighted sum 
of the spectral distribution for the cr band and for the tt 
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band, each of the two calculated according to Eq. ij^J: 

Ntot{r,T) = w^N„{f,E)+w^N^{r,E) (36) 

For the special case of MgB2, the weight of the spec- 
tral distribution of the tt band is given by Wtt = 
7V(0)(^V(^(0)^'"^ + iV(0)('^)) = 0.577, which implies 

= 1 - = 0.423 (e.g. see 
In Fig. El we show numerical results for the two band 
quasiparticle spectrum in the special case of MgB2. The 
total quasiparticle spectrum exhibits much more struc- 
ture due to the superposition of the spectra of the two 
single bands. 



B. The analytical Gap Model for Two Bands 

In order to extend the analytical gap model from chap- 
ter ljIIB|l to the case of a two gap superconductor, one has 
to introduce two separate gap models for the pairing po- 
tential, one for each band, and each of the two in the form 
of Eq. 10. Hence, one has to introduce two lengthscales, 

one for each of the two bands: ^i"^^ and . Following 
the discussion in section ijlll A|l . the calculations carried 
out in chapter ljIIB|l have then to be repeated for each of 
the two bands separately in order to obtain the quasiclas- 
sical Green's function for both bands: /^(r, kp,iSn) and 
fTr{r,kp,iSn)- A coupling of the two bands first occurs 
in the gap equation Ij34|l . 

In full analogy to the procedure in chapter ljIIB|l . the 
Green's function /(r, kp, *£n) can be calculated for each 
of the two bands. The two band gap equation l(3ijl can 
again be differentiated at the center of the vortex xq = 
with respect to the position a;o and determines the two 
lengthscales S^i"^^ and <^i'^^ for the full temperature range: 

= 2TrT^Xaa' ^ {dxofa'ixQ =0))FS.a' 
i'u a' 0<e„<iJc 

(37) 

In Fig. we show the results from these calculations, 
again compared to the results from a fully self-consistent 
solution of the two band gap equation. As in the single 
gap case, the linear profile for the pairing potential re- 
sults in a qualitatively correct behaviour of ^i'^'' (T) and 
^i'^\T) for the full temperature range. The lengthscales 
and ^i^^ diverge for T ^ Tc and vanish for T ^ 0. 
Thus, we find a Kramer-Pesch effect in both bands of a 
clean two band system. Again, for low temperatures the 
linear profile matches nicely the self-consistently calcu- 
lated curve but shows deviations for higher temperatures. 
As in the single band case, the lengthscales S^u''^ and S^v'^ 
bend towards ^i"^ — S^v'^ = not until very low tem- 
peratures 0.05 Tc) and thus one should again be care- 
ful with extrapolations from higher temperatures towards 
T = (see inset in Fig-CJ. 




E/keTc 



Figure 6: The quasiparticle spectrum in the two band case, 
again for several distances from the vortex core, (a) shows the 
quasiparticle spectrum of the cr band and (b) that of the tt 
band. In (c), the total quasiparticle spectrum is being plotted. 
As before, the distances span from r = (rearmost curve) to 
r = 5 ^oo (foremost curve) in steps of 0.2 and an imaginary 
part of 5 = 0.1 was used. 
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Figure 7: The lengthscales (,v (T) and 5„ (T) obtained from 
a linear profile of the pairing potential (solid line) and from 
a numerical self-consistent solution of the gap equation (di- 
amonds with dashed guidelines). Inset: Blow-up of the low 
temperature region. 
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Figure 8: The numerical results for the vortex core radii 
as a function of temperature for the a band (diamonds) and 
the TT band (stars) calculated in the dirty limit with a ratio of 
the diffusion constants V^"^ = Q.2V^''\ One can clearly see a 
saturation for low temperatures. The lines are guides for the 
eye. 

We also show the comparison with a dirty limit cal- 
culation of the two band gap equation using the Usadel 
formalism introduced in section IlII see Fig.|Hl In this 
case we used slightly different parameters to calculate the 
coupling matrix Xaa' than in the clean limit calculations. 
For the ratio of the diffusion constants we have chosen 
the value proposed by Koshelev and Golubov in as 
pC*^) = 0.21)^^^ leading to a ratio of the coherence lengths 
to be = 0A77^!^\ As expected, we find a finite core 
size for both bands even for low temperatures and thus, 
as in the single gap case, no Kramer-Pesch effect exists 
in the dirty limit. 

In order to compare the behaviour of the two bands, it 
is useful to study the ratio ^i^^ /^i'^^ In Fig. El we show 
this ratio for the clean and the dirty limit. The clean limit 
results were obtained, on the one hand, from a numeri- 
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Figure 9: The ratio £,v /(,v for the clean and the dirty limit. 
Solid line: results from the linear profile of the pairing poten- 
tial for the clean limit. Diamonds with dashed guideline: nu- 
merical solution of the gap equation in the clean limit. Stars 
with dotted guideline: Usadel formalism results for the dirty 
limit. 

cal solution of the gap equation and, on the other hand, 
from the analytical model for the pairing potential. The 
dirty limit results stem from calculations within the Us- 
adel formalism. For high temperatures (T Tc), both 

^■i'^'' and ^i'^'' diverge and the ratio ^i^V^"'^^ tends to 1 
in the clean as well as in the dirty limit. For low tem- 
peratures (T ^ 0), the clean and the dirty limit results 
differ greatly. In the clean limit, both lengthscales van- 
ish, however, in a way that leads to a strong growth of 
the ratio ^i'^^ i.e. the vortex core size in the tt band 

decreases more slowly as a function of temperature than 
the one in the a band. The agreement between the nu- 
merical results and those from the linear profile is very 

(tt) 

good here. In the dirty limit, the saturation of both 
and leads to a saturation of the ratio ^i^^ /S^v'' ■ 

The quasiparticle spectrum for two dirty bands is a 
simple superposition of two single band spectra as can 
be seen in Fig. [Till 

IV. INDUCED KRAMER-PESCH EFFECT IN A 
DIRTY TT BAND 

In present high-quality samples of MgB2 it is believed 
that the intraband scattering rate in the a bands is 
smaller than the gap value, while in the tt bands the 
intraband scattering rate is larger than the correspond- 
ing gap value ^21. .22. . Therefore we want to discuss 
in this section - as first suggested by M. Eschrig j23| - a 
model with a clean a band and a dirty tt band. Since the 
gap profiles in the two bands can not vary independently 
as they are coupled via the gap equation, we expect that 
for low temperatures they show the same asymptotic be- 
haviour. But it is not at all clear if we find a decrease of 
and for T ^ as in the clean limit calculations 
or a saturation value as in the dirty limit approach. To 
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Figure 10: The quasiparticle spectrum in the two band dirty 
limit approach, again for several distances from the vortex 
core. As before, the distances span from r = (rearmost 
curve) to r = 10 ^cso (foremost curve) in steps of O.S^oo. 




Figure 11: The numerical results for the votex core radii 
as a function of temperature for the clean a band (diamonds) 
and the dirty n band (stars) calculated in a mixed model 
with equal coherence lengths in the two bands = 
The clean a band and the dirty tt band show both a signifi- 
cant Kramer-Pesch effect for low temperatures. The ratio of 
the core radii in the tt and in the a band as a function of 
temperature calculated for a clean a and a dirty tt band is 
shown in the inset. Again in both plots the lines are guides 
for the eye. 




obtain a result for the two vortex core radii we have to 
solve the two band gap equation where we calcu- 
late fa from the Riccati amplitudes and ba according 
to while the Fermi surface averaged Green's func- 
tion (/7r)_FS,ir is replaced by the Usadel Green's function 
i^.n.(r, £„) that is a solution of a simple boundary value 
problem as described in section IlII C|l . Since we have two 
different lengthscales in the clean and the dirty band case 



Figure 12: The quasiparticle spectrum in the two band case, 
again for several distances from the vortex core. Plot (a) 
shows the quasiparticle spectrum of the clean a band and 
(b) that of the dirty tt band. In (c), the total quasipar- 
ticle spectrum is being plotted for the mixed model. As 
before, the distances span from r = (rearmost curve) to 
r = 10 ^oo (foremost curve) in steps of ^oo and an imaginary 
part of 5 = 0.1 was used. 
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Coo = -^xlj and = Y ^j^^ we have to fix the ratio of 

the two coherence lengths which in general will depend 
on the sample quality, i.e. the scattering rate in the tt 
band. Assuming the diffusion constant P^^-' to be of the 
order of 

we chose the ratio to be ^^V'foo'' — 1- We do not expect 
the essential qualitative results to depend on the special 
value of this ratio which we have checked numerically for 
a different ratio. In Fig. we show the numerical results 
of this calculation. In our calculation we find a distinct 
decrease of the vortex core size for low temperatures even 
in the dirty tt band, that can be understood as an induced 
Kramer-Pesch effect due to the real space coupling of the 
two gap amplitudes via the gap equation. The infiuence 
of the dirty tt band is refiected in a larger ratio of ^^.^■'/Ci''^' 
for low temperatures than in the pure clean limit or the 
pure dirty limit approach. 

In Fig. El we show the quasiparticle spectrum for the 
mixed model, consisting of a clean a band and a dirty tt 
band. Looking at the a and tt band spectra individually, 
it is obvious that the zero energy bound state only exists 
in the clean band. Thus, a zero energy peak is visible 
in the total quasiparticle spectrum, but it is lifted by an 
offset due to the fiat spectrum of the dirty band. Re- 
stricting to the vortex core region, this leads to a picture 
that is very similar to the clean single band case. 

V. CONCLUSIONS 

We have studied the Kramer-Pesch effect in both a sin- 
gle gap and a two gap superconductor using parameters 
relevant for MgB2. We have presented a model for the 



vortex core that allows an analytical solution for the Ric- 
cati equations in both cases and compared it with fully 
self-consistent numerical solutions. This model is useful 
for approximate calculations, in particular at low tem- 
peratures. 

In the two gap case we find that the Kramer-Pesch 
effect is present in both bands in the clean limit. At 
high temperatures the sizes of the two components of 
the vortex core become equal, while at low temperatures 
the size of the vortex core in the tt band with the smaller 
gap decreases more slowly with decreasing temperature. 

We also investigated the dirty limit within the frame- 
work of the Usadel equations. In a single band supercon- 
ductor there is no Kramer-Pesch effect in the dirty limit 
as has been discussed before Q. However, in a two gap 
superconductor an interesting new situation arises, when 
one band is in the clean limit while the other is in the 
dirty Hmit, as is believed to be the case in high quality 
MgB2 samples. In this particular case we find that even 
in the dirty band a Kramer-Pesch effect exists, induced 
by the clean band. 

Traditionally the Kramer-Pesch effect is difficult to ob- 
serve experimentally due to the disturbing effect of im- 
purities. Our calculations indicate that there is a better 
chance of observing this effect in MgB2 by scanning tun- 
neling microscopy (STM) imaging of the vortex core as a 
function of temperature, because the tt band may be in 
the dirty limit and still showing a decreasing vortex core 
size with decreasing temperature, as seen in Fig. ITTl 
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